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Abstract 


's^ = 2.76 TeV was studied in the 


Direct photon production at mid-rapidity in Pb-Pb collisions at 
transverse momentum range 0.9 < px < 14 GeV /c. Photons were detected with the highly segmented 
electromagnetic calorimeter PHOS and via conversions in the ALICE detector material with the 
pair reconstructed in the central tracking system. The results of the two methods were combined and 
direct photon spectra were measured for the 0-20%, 20^0%, and 40-80% centrality classes. For all 
three classes, agreement was found with perturbative QCD calculations for pj ^5 GeV/c. Direct 
photon spectra down to px ~ 1 GeV /c could be extracted for the 20-40% and 0-20% centrality 
classes. The significance of the direct photon signal for 0.9 < px < 2.1 GeV /c is 2.6(7 for the 0- 
20% class. The spectrum in this px range and centrality class can be described by an exponential 
with an inverse slope parameter of (297 ± 12*“±4P^'‘*)MeV. State-of-the-art models for photon 
production in heavy-ion collisions agree with the data within uncertainties. 
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1 Introduction 


The theory of the strong interaction, Quantum ChromoDynamics (QCD), predicts a transition from or¬ 
dinary nuclear matter to a state where quarks and gluons are no longer confined to hadrons |[ll0!- The 
creation and study of this deconfined partonic state, the Quark-Gluon Plasma (QGP), is the major objec¬ 
tive in the experimental program of heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) 
M and the Large Hadron Collider (LHC) iffl]. 

Direct photons, defined as photons not originating from hadron decays, are a valuable tool to study de¬ 
tails of the evolution of the medium created in heavy-ion collisions. Unlike hadrons, direct photons 
are produced at all stages of the collision and escape from the hot nuclear matter basically unaffected 
I lol . delivering direct information on the conditions at the time of production: prompt direct photons 
produced in hard scatterings of incoming partons provide information on parton distributions in nuclei; 
deconfined quark-gluon matter as well as hadronic matter created in the course of the collision emit ther¬ 
mal direct photons, carrying information about the temperature, collective flow and space-time evolution 
of the medium il7n . Different transverse momentum (px) regions are dominated by photons emitted 
at different stages of the collision. Prompt direct photons follow a power law spectrum and dominate 
at high transverse momentum (px ^ 5 GeV/c). At lower transverse momenta (px < 4 GeV/c) one 
expects contributions from the thermalized partonic and hadronic phases with an approximately expo¬ 
nential spectrum ifl^ H^ . In addition, other direct photon production mechanisms, like the interaction 


of hard scattered partons with the medium (”jet-photon conversion”) Il2d. 12111 . may be important for 
Px < lOGeV/c. 

The direct photon spectrum at low px, therefore, contains information on the initial temperature and 
space-time evolution of the thermalized medium created in heavy-ion collisions. The observed thermal 
direct photon spectrum is a sum of contributions from all stages of the collision after thermalization, 
where the earliest, hottest stage and later, cooler stages can make comparable contributions 0. High 
photon emission rates at the largest temperatures in the early stage are compensated by an expanded 
space-time volume and blue-shift due to radial flow in the later stage. This complicates the interpretation 
of inverse slope parameters of direct photon spectra, but a correlation between the slope and the initial 
temperature still exists Q. 

The first measurement of a direct photon spectrum in relativistic A-A collisions was presented by the 
WA98 collaboration 11241] . The direct photon yield was measured at the CERN SPS in central Pb-Pb 
collisions at =17.3 GeV in the range 1.5 < px < 4 GeV/c. The signal can be interpreted either 
as thermal photon radiation from a quark-gluon plasma and hadronic gas or as the effect of multiple 
soft scatterings of the incoming partons without the formation of a QGP The PHENIX experiment 
measured the direct photon spectrum in Au-Au collisions at = 200 GeV in the range 1 < px ^ 


20 GeV/c i25l. 12611. It was found that at high px (5 < px ^ 21 GeV/c) the direct photon spectrum 


measured in Au-Au collisions agrees with the one measured in pp collisions at the same energy after 
scaling with the number of binary nucleon-nucleon collisions (/7coii)- Scaling of high-px direct photon 
production with Vcoii in Pb-Pb collisions at LHC energy was confirmed by the ATLAS 1271] and CMS 
ll28n experiments in the measurement of isolated photons, i.e., photons with little hadronic energy in a 
cone around them, in the ranges 22 < px < 280 GeV/c and 20 < px < 80 GeV/c, respectively. The 
absence of suppression of high pj isolated photons in A-A collisions with respect to Vcoii scaled pp 
collisions, in contrast to the observed suppression of hadrons, is consistent with the latter being due to 
energy loss of hard scattered quarks and gluons in the medium. 

Direct photon production at low px (< 3 GeV/c) in Au-Au collisions at = 200 GeV was studied by 
the PHENIX experiment in the measurement of virtual photons {e^e^ pairs from internal conversions) 


12911 and with real photons t30ll . A clear excess of direct photons above the expectation from scaled 
pp collisions was observed. The excess was parameterized by an exponential function with inverse 
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slope parameters Teff = 221 ± 19®*^*^ ± 19®^®^ MeV (virtual photon method 1291] ) and Ten = 239 ± 25 


:stat. 


V^yst MeV (real photon method 1301]) for the 0-20% most central collisions. The measured spectrum 
can be described by models assuming thermal photon emission from hydrodynamically expanding hot 
matter with initial temperatures in the range 300 — 600 MeV Bin . The measurement of a direct-photon 
azimuthal anisotropy (elliptic flow), which was found to be similar in magnitude to the pion elliptic 
flow at low pt in Au-Au collisions at = 200 GeV 11321] . provides a further important constraint for 
models. The simultaneous description of the spectra and elliptic flow of direct photons currently poses a 
challenge for hydrodynamic models Q. 

In this letter, the first measurement of direct photon production for < 14 GeV/c in Pb-Pb collisions 
at 


A]]]];]' = 2.76 TeV is presented. 


2 Detector setup 

Photons were measured using two independent methods: by the Photon Conversion Method (PCM) 
and with the electromagnetic calorimeter PROS. In the conversion method, the electron and positron 
tracks from a photon conversion were measured with the Inner Tracking System (ITS) and/or the Time 
Projection Chamber (TPC). 


The ITS 13411 consists of two layers of Silicon Pixel Detectors (SPD) positioned at a radial distance of 
3.9 cm and 7.6 cm, two layers of Silicon Drift Detectors (SDD) at 15.0 cm and 23.9 cm, and two layers of 
Silicon Strip Detectors (SSD) at 38.0 cm and 43.0 cm. The two innermost layers cover a pseudorapidity 
range of |t]| <2 and |t]| < 1.4, respectively. The TPC i35n is a large (85 m^) cylindrical drift detector 
filled with a Ne-C02-N2 (90-10-5) gas mixture. It covers the pseudorapidity range |t] | < 0.9 over the full 
azimuthal angle with a maximum track length of 159 reconstructed space points. With the magnetic field 
of B = 0.5 T, and tracks can be reconstructed down to px ~ 50 MeV/c, depending on the position 
of the conversion point. The TPC provides particle identification via the measurement of the specific 
energy loss (dE/dx) with a resolution of 5.2% in pp collisions and 6.5% in central Pb-Pb collisions S]. 
The ITS and the TPC were aligned with respect to each other to the level of less than 100 pm using 
cosmic-ray and pp collision data OTP . Particle identification is furthermore provided by the Time-of- 
Flight (TOE) detector i38n located at a radial distance of 370 < r < 399 cm. This detector consists of 
Multigap Resistive Plate Chambers (MRPC) and provides timing information with an intrinsic resolution 
of 50 ps. 


PROS i39h is an electromagnetic calorimeter which consists of three modules installed at a distance of 
4.6 m from the interaction point. It subtends 260° < (p < 320° in azimuth and |t]| < 0.13 in pseudo¬ 
rapidity. Each module consists of 3584 detector cells arranged in a matrix of 64 x 56 lead tungstate 
crystals each of size 2.2 x 2.2 x 18 cm^. The signal from each cell is measured by an avalanche photodi¬ 
ode (APD) associated with a low-noise charge-sensitive preamplifier. To increase the light yield, reduce 
electronic noise, and improve energy resolution, the crystals, APDs, and preamplifiers are cooled to a 
temperature of —25 °C. The resulting energy resolution is Oe/E = {\.3%/E) © (3.3%/\/E) © 1.12%, 
where E is in GeV. The PROS channels were calibrated in pp collisions by aligning the 71^ peak position 
in the two-photon invariant mass distribution. 


Two scintillator hodoscopes (VO-A and VO-C) t40ll subtending 2.8 < T] < 5.1 and —3.7 < T] < —1.7, 
respectively, were used in the minimum bias trigger in the Pb-Pb run. The sum of the amplitudes of 
VO-A and VO-C served as a measure of centrality in the Pb-Pb collisions. 


3 Data analysis 

This analysis is based on data recorded by the AEICE experiment in the first ERG heavy-ion run in the 
fall of 2010. The detector readout was triggered by the minimum bias interaction trigger based on trigger 
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signals from the VO-A, VO-C, and SPD detectors. The efficiency for triggering on a Pb-Pb hadronic 
interaction ranged between 98.4% and 99.7%, depending on the minimum bias trigger configuration. The 
events were divided into centrality classes according to the VO-A and VO-C summed amplitudes. Only 
events in the centrality range 0-80% were used in this analysis. To ensure a uniform track acceptance 
in pseudorapidity rj, only events with a primary vertex within ±10 cm from the nominal interaction 
point along the beam line (z-direction) were used. After offline event selection, 13.6 x 10^ events were 
available for the PCM analysis and 17.7 x 10^ events for the PHOS analysis. 


The direct photon yield is extracted on a statistical basis from the inclusive photon spectrum by com¬ 
paring the measured photon spectrum to the spectrum of photons from hadron decays. The yield of 
7r*’s, which contribute about 80 — 85% of the decay photons (cf. Fig. [T]l, was measured simultaneously 
with the inclusive photon yield. Besides photons from 71° decays, the second and third most important 
contributions to the decay photon spectrum come from rj and ft) decays. 


An excess of direct photons above the decay photon spectrum can be quantified by the pj dependent 
double ratio 

Tmcl / /decay _ /ncl 


Ry — 




0 

pai'am 




0 


param 


V ■ 

/decay 

where /inci is the measured inclusive photon spectrum, TTp^ram ^ parameterization of the measured 71° 
spectrum, andj^ecay the calculated decay photon spectrum. The PCM and PHOS 71 ° measurements are 
described in 14IH . The double ratio has the advantage that some of the largest systematic uncertainties 
cancel partially or completely. Using the double ratio, the direct photon yield can be calculated from the 
inclusive photon yield as 
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/direct — /incl /decay — (1 d ) ' ^"cl 
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The PCM and PHOS analyses were performed independently. Combined direct photon spectra were 
determined based on combined double ratios and combined inclusive photon spectra. In contrast to 
taking the average of the PCM and PHOS direct-photon spectra, this approach allowed us to use the 
information from both measurements also when one measurement of Ry fluctuated below unity. 

In the PCM analysis, photons are reconstructed via a secondary vertex finding algorithm which provides 
displaced vertices with two opposite-charge daughters. The positively and negatively charged daughter 
tracks are required to contain reconstructed clusters in the TPC. Only tracks with a transverse momentum 
above 50 MeV/c and a ratio of the number of reconstructed TPC clusters over the number of findable TPC 
clusters (accounting for track length, spatial location and momentum) larger than 0.6 were considered. 
To identify and e^, the specific energy loss in the TPC 0611 was required to be within a band of 
[—3a,5a] around the average electron dE/dx, and be more than 3a above the average pion dE/dx 
(where the second condition is only applied for tracks with pj > 0.4 GeV/c). Tracks with an associated 
signal in the TOE detector were only accepted as electron candidates if they were consistent with the 
electron hypothesis within a ±5 a band. The vertex finding algorithm uses the Kalman filter technique 
for the decay/conversion point and four momentum determination of the neutral parent particle (U°) 14211 . 
U°s result from 7 conversions but also from strange particle decays (K°, A or A). Further selection was 
performed on the level of the reconstructed U°. U°s with a decay point with radius r < 5 cm were rejected 
to remove 71 ° and rj Dalitz decays. The transverse momentum component qj = Pe sin dvo,e H of the 
electron momentum, pe, with respect to the U° momentum was restricted to qj < 0.05 GeV/c. Based 
on the invariant mass of the pair and the pointing of the U° to the primary vertex, the vertex finder 
calculates a X^(/) value which reflects the level of consistency with the hypothesis that the U° comes 
from a photon originating from the primary vertex. A selection based on this X^{y) value was used to 
further reduce contamination in the photon sample. Random associations of electrons and positrons were 
further reduced by making use of the small opening angle of the pair from photon conversions at 
the conversion point. 
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The raw photon spectrum, constructed from the secondary vertex candidates passing the selection de¬ 
scribed above, was corrected for the reconstruction efficiency, the acceptance and the contamination. The 
detector response was simulated for Pb-Pb collisions using HIJING 11441] together with the GRANT 3.21 
transport code 1451] . The resulting efficiency correction is dominafed by fhe conversion probabilify of 
phofons in fhe ALICE material. The integrated maferial budge! of fhe beam pipe, fhe ITS and fhe TPC 
for r < 1.8 m corresponds fo (11.4 ±0.5)% of a radiation lengfh Aq, resulting in a photon conversion 
probabilify thaf safurafes af abouf 8.5% for Pt:>2 GeV/c IStiL l42h . The photon finding efficiency for 
converted phofons is of fhe order of 50-65% over fhe measured pj range for all cenfralifies. The purify 
of fhe photon candidafe sample for pt < 3 GeV/c exfracfed from simulation is 98-99% in peripheral 
and 91-97% in the most central collisions. Eurthermore, secondary photon candidates, mainly photons 
from the decay —>• 271** —>• 4 7 , not removed by the selection, were subtracted statistically based 
on the measured spectrum 14611 . A correction of less than 2% for photons from pile-up collisions was 
applied for the 40-80% class for pj <2 GeV/c. At higher pj and for more central classes this correction 
is negligible. 


Centrality 

0 -20% 

20 =10% 

40- 

-80% 

px (GeV/c) 

1.2 

5.0 

1.2 

5.0 

1.2 

5.0 

7 u„;i yield 

Track quality (A) 

0.6 

0.6 

0.2 

0.2 

0.2 

0.7 

Electron PID (A,B) 

1.5 

6.9 

0.9 

4.8 

0.7 

4.0 

Photon selection (A,B) 

4.0 

1.8 

2.4 

2.1 

1.5 

1.3 

Material (C) 

4.5 

4.5 

4.5 

4.5 

4.5 

4.5 

YtDCl/jt^ 

Track quality (A) 

0.7 

1.7 

0.8 

0.4 

0.6 

1.3 

Electron PID (A,B) 

1.2 

4.8 

0.9 

3.8 

0.9 

4.0 

Photon selection (A,B) 

3.2 

3.2 

3.0 

1.5 

2.5 

2.4 

71*’ yield (A) 

1.6 

2.9 

1.7 

2.7 

0.5 

3.0 

Material (C) 

4.5 

4.5 

4.5 

4.5 

4.5 

4.5 

Ydecay /l^ 

71*’ spectrum (B) 

0.5 

1.2 

0.8 

1.8 

0.5 

3.2 

T] yield (C) 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

T] shape (B) 

1.6 

0.5 

1.2 

0.2 

1.0 

0.2 

Total Ry 

6.2 

8.1 

5.7 

7.0 

5.7 

8.3 

Total 7inci 

6.2 

8.5 

5.2 

6.9 

4.8 

6.2 


Table 1: Summary of the systematic uncertainties of the PCM analysis in percentage. Uncertainties are character¬ 
ized according to three categories: point-by-point uncorrelated (A), correlated in pj with magnitude of the relative 
uncertainty varying point-by-point (B), and constant fractional uncertainty (C). Items in the table with categories 
(A,B) summarize sources of uncertainties which are either of type A or B. 


In the PHOS analysis, clusters (each cell of the cluster must have at least one common edge with another 
cell of the cluster) were used as photon candidates. To estimate the photon energy, the energies of cells 
with centers within a radius Rcoie = 3.5 cm from the cluster center of gravity were summed. Compared 
to the full cluster energy, this core energy (Ecore) is less sensitive to overlaps with low-energy clusters in a 
high multiplicity environment. The non-linearity in the conversion of the reconstructed to the true photon 
energy introduced by this approach is reproduced by GEANT3 Monte Carlo simulations. The contribu¬ 
tion of hadronic clusters was reduced by requiring E^ter >0.3 GeV, Vceiis > 2 and by accepting only 
clusters above a minimum lateral cluster dispersion 14111 . The latter selection rejects hadrons punching 
though the crystal and producing a large signal in the photodiode of a single cell. With a minimum time 
between bunch crossings of 525 ns, possible pile-up contributions from other bunch crossings is removed 
by a loose cut on the cluster arrival time |f| < 150 ns. Eor systematic uncertainty studies, photons were 
also reconstructed with a px-dependent dispersion cut and with a charged particle veto (CPV) cut on the 
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Centrality 

0 - 20 % 

20 ^ 0 % 

40- 

-80% 

Pj (GeV/c) 

2 

10 

2 

10 

2 

10 

rinci yield 

Efficiency (B) 

3.0 

3.0 

0.7 

0.7 

2.5 

2.5 

Contamination (B) 

2.0 

2.0 

1.3 

1.3 

2.9 

0.5 

Conversion (C) 

1.7 

1.7 

1.7 

1.7 

1.7 

1.7 

Acceptance (C) 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

^Global E scale (B) 

9.6 

9.0 

6.1 

5.9 

5.8 

6.3 

*Non-linearity (B) 

2.2 

0.1 

2.1 

0.1 

2.0 

0.1 

TT® yield 

Yield extraction (A) 

2.7 

4.0 

3.1 

5.2 

1.8 

2.9 

Efficiency (B) 

1.8 

1.8 

2.7 

2.2 

2.5 

2.5 

Acceptance (C) 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Pileup (C) 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Feed-down (B) 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

T^ecay/ ^ 

71*’ spectrum (B) 

1.3 

4.3 

1.8 

1.8 

1.8 

1.8 

T] contribution (B) 

2.2 

1.7 

2.2 

1.6 

2.1 

1.6 

Total Ry 

6.8 

7.9 

5.9 

6.5 

6.1 

6.0 

Total 7nci 

12.4 

12.7 

9.7 

10.0 

9.8 

9.6 


Table 2: Summary of systematic uncertainties of the PHOS analysis in percentage. Uncertainties are characterized 
according to three categories: point-by-point uncorrelated (A), correlated in pj with magnitude of the relative 
uncertainty varying point-by-point (B), and constant fractional uncertainty (C). Uncertainties marked with * cancel 
in the double ratio Ry. 


distance between the PHOS cluster position and the TOsition of extrapolated charged tracks on the PHOS 
surface to suppress clusters from charged particles 14 ill . Both dispersion and CPV cuts were tuned using 
pp collision data to provide a photon efficiency at the level of 96-99%. 

The product of acceptance and efficiency (A • e) was estimated by embedding simulated photon clus¬ 
ters into real events and applying the standard reconstruction. PHOS properties (energy and position 
resolutions, residual de-calibration, absolute calibration, non-linear energy response) were tuned in the 
simulation to reproduce the pj dependence of the tt** peak position and width 14111 . In peripheral events, 
A • £ for the default selection (no dispersion cut, no CPV cut) has a value of about 0.022 at px = 1 GeV/c. 
For higher pT^ A • £ decreases and saturates at about 0.018 for px ^ 5 GeV/c. The decrease of A • £ with 
Pj results from the use of Ecore- In central collisions, A • £ increases by up to about 10% due to clus¬ 
ter overlaps. Applying the dispersion and CPV cuts, the efficiency is reduced by 5-10% in peripheral 
collisions and the centrality dependence becomes negligible. 

The contamination of the photon spectrum measured with PHOS originates mainly from and p, h an¬ 
nihilation in PHOS, with other contributions being much smaller. Application of the dispersion and CPV 
cuts reduces the overall contamination at px ~ 1-5 GeV/c from about 15% to 2-3% and down to 1-2% at 
Pj ~ 3^ GeV/c. The subtraction of contamination is based on a data driven approach: the probability to 
pass the CPV and dispersion cuts and the calorimeter response to hadrons are estimated using identified 
71*, p tracks; the photon candidate spectra, measured with different cuts (default, dispersion, CPV, both) 
were decomposed into 7 , tt*, p and h contributions, assuming equal contamination from p and h. The 
contamination calculated in this way agrees with that estimated from a HHING simulation. Finally, the 


photon contribution from K\ 
1461 as in the PCM analysis. 


271'^ —)■ 47 decays was subtracted based on the measured /ff spectrum 


To calculate the 7decayratio, a Monte Carlo approach was used to simulate particle decays into photons 
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both for the PCM and the PROS analysis. The largest contributions come from Ti^, rj, and ft) decays. 
Contributions of other hadrons were also included but were found to be negligible. To allow for a 
cancellation of some uncertainties common to the photon and yield in Eq. ([T), each analysis (PCM, 
PROS) used the spectrum measured with the respective method. 



10 

(GeV/c) 


Fig. 1: (Color online) Relative contributions of different hadrons to the total decay photon spectrum as a function 
of the decay photon transverse momentum (PCM case). 


The T] meson contribution is estimated by using two approaches which assume: (i) transverse mass (mj) 
scaling of the 71** and the rj spectrum which is consistent with measurements at RRIC SStI] or (ii) that 
the pt spectrum of the rj has the same shape as the spectrum S as both particles should be affected 
by radial flow in the same way due to their similar masses. The maximum deviation between these two 
cases occurs at px ~ 2.5 GeV/c where (i) corresponds to a 77 / 71 ° ratio of about 0.4 whereas (ii) gives a 
ratio of about 0.5. The absolute yield of rj mesons in both cases was fixed at px > 5 GeV/c to reproduce 
the measured 77/71° ratio at = 200 GeV: 0.46 ± 0.05 14811 . The statistical precision of the 77 signal 


in the 2010 and 2011 data sets is too low to further constrain these two assumptions with a measurement 
of the 77 spectrum. The average of these two cases is used for the decay photon calculation, while half 
the difference is taken as a contribution to the systematic uncertainty of the 77 meson contribution in 
addition to the normalization uncertainty quoted above. The contribution of co meson decay photons is 
below ~ 3% and mj scaling of the measured 71° spectrum with (dAcB/d/ 77 x)/(dA;jo/d/ 77 x) = 0.9 is used 
ll^ . The relative contributions of the different hadrons to the total decay photon spectrum are shown in 
Eig.ffl 


The main sources of systematic uncertainties in the determination of the inclusive photon spectrum and 
Ry for the PCM analysis are listed in Table [T] The two largest uncertainties are related to the material 
budget of the AEICE detector and the Monte Carlo-based efficiency corrections to the inclusive pho¬ 
ton and 71° spectra. The material budget uncertainty was estimated in pp collisions by comparing the 
measured number of converted photons (normalized to the measured charged particle multiplicity) with 
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GEANT simulation results in which particle yields from PYTHIA and PHOJET were used as input. Un¬ 
certainties related to track selection and electron identification were estimated by variation of the cuts. 
Eor instance, we observe a small variation in results depending on the minimum threshold for electron 
tracks. This is most likely related to different tracking performance for real data and in the Monte Carlo 
simulation for low-pT particles (px ^ 50 MeV/c). The uncertainty related to the choice of this threshold 
was estimated by increasing the minimum px from 50 MeV/c up to 100 MeV/c. Uncertainties related 
to falsely reconstructed electron-positron pairs from Dalitz decays as conversion pairs were obtained by 
varying the minimum radial distance of reconstructed electron tracks from the standard value of 
^min = 5 cm up to /?min = 10 cm. The estimation of the systematic uncertainty of the electron selection 
includes a contribution estimated by the variation of the dE/dr cuts. 


In the double ratio Ry, many uncertainties partially cancel. The uncertainties on Ry were therefore ob¬ 
tained by evaluating the effect of cut variations directly on Ry. Uncertainties related to the decay photon 
spectrum are similar for PCM and PHOS analyses: they include the uncertainty due to the 7i^ spectrum 
parameterization, difference of shapes of spectra measured by PCM and PHOS, and uncertainties due 
to the shape and absolute normalization of the Tj spectrum. Uncertainties due to contributions of other 
hadrons are negligible. 


The main systematic uncertainties of the PHOS analysis are summarized in Table |2] Eor the inclusive 
photon spectrum, the uncertainty of the efficiency calculation is estimated comparing the PID cut effi¬ 
ciency in Monte Carlo and real data. The contamination uncertainty is estimated comparing the photon 
purity calculated with a data driven approach and with Monte Carlo HHING simulations. The conversion 
probability is estimated comparing yields in pp collisions with and without magnetic field. The global 
energy and non-linearity uncertainties, which mostly cancel in Ry, are estimated comparing calibrations 
based on the tt** peak position and on the electron E /p peak position. The centrality dependence of 
the energy scale uncertainty results from the larger background under the 71^ peak in central events and 
therefore larger uncertainties in the peak position. 


A more detailed description of the single photon selection and e^cially of the additional uncertain¬ 
ties for both the PCM and PHOS analyses can be found in Ref. 84 111 . 


The comparison of the individual PHOS and PCM inclusive photon spectra, normalized to the averaged 
spectrum, is shown in Pig. |2] Statistical and point-to-point uncorrelated systematic uncertainties (type A) 
are combined and presented as error bars, point-to-point correlated systematic uncertainties (type B) are 
shown as boxes, and common normalization systematic uncertainties (type C) are shown as bands around 
unity. The uncertainties are dominated by px-correlated contributions. The individual PHOS and PCM 
double ratios are shown in Pig. |3] The partial cancellation of the energy scale uncertainties (PHOS) and 
the material budget uncertainties (PCM) is taken into account in the shown uncertainties. 


The level of agreement between the PHOS and PCM inclusive photon spectra and double ratios Ry was 
quantified taking into account the correlation of the uncertainties in pj and centrality. To this end, pseudo 
data points for the ratio of the PHOS and PCM inclusive photon spectra and double ratios were generated 
simultaneously for all three centrality classes under the assumption of the null hypothesis that the ratio 
is unity for all points, i.e., that both measurements result from the same original distribution. The type B 
and C systematic uncertainties give rise to a shifted baseline, around which the pseudo data points are 
drawn from a Gaussian with a standard deviation given by the statistical and type A uncertainties. A 
test statistic t was defined as the sum of the squared differences of the pseudo data points with respect 
to the null hypothesis in units of the type A and statistical uncertainties. A p-value was calculated 
as the fraction of pseudo experiments with values of t larger than observed in the real data 115011 . The 
corresponding significance in units of the standard deviation of a one-dimensional normal distribution 
was calculated based on a two-tailed test. The PHOS and PCM inclusive photon spectra were found 
to agree within 1.2 standard deviations, the PHOS and PCM double ratios agree within 0.4 standard 
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Fig. 2: (Color online) Comparison of inclusive photon spectra measured with PCM and PHOS in the 0-20%, 20- 
40%, and 40-80% centrality classes. The individual spectra were divided by the corresponding combined PCM 
and PHOS spectrum. The shown errors only reflect the uncertainties of the individual measurements. The boxes 
around unity indicate normalization uncertainties (type C). 


deviations. 


4 Results 

The inclusive photon spectra and double ratios of the PCM and PHOS analyses are combined as two 
independent measurements to obtain the error-weighted average. The uncertainties common to both 
measurements (trigger efficiency, centrality determination, etc.) are negligible in comparison to the 
uncorrelated, analysis-specific uncertainties. For each cenfralify selection fhe average double rafio Ry is 
used fogefher wifh fhe averaged inclusive photon speefrum fo obfain fhe final direef photon speefrum, 
according to Eq. |2j For fhe 0-20% cenfralify class and pt = 2 GeV/c, fhis resulfs in fype A, B, and C 
sysfemafic uncerfainfies of a a = 2.5%, rte = 2.3%, and Gq = 3.0% for fhe combined double ratio and 
of a A = 20%, rtB = 18%, rtc = 24% for fhe combined direef photon speefrum. 

The combined PCM and PHOS double ratios Ry measured for fhree cenfralify classes are shown in 
Fig. m A direef photon excess is observed for all cenfralify classes for pj ^4 GeV/c, and also for 
1 ~ Pt ^ 4 GeV/c in fhe mosf cenfral class. The measuremenfs are compared wifh fhe expeefed Ry 
for fhe prompf photon confribufion as calculafed wifh nexf-fo-leading-order (NFO) perturbafive QCD 
calculations. The prompf photon expeefafions in Fig. |4] were defermined as 1 +A^co1i7pQCd// decay where 
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Fig. 3: (Color online) Comparison of double ratios Ry measured with PCM and PHOS for the 0-20%, 20^0%, and 
40-80% centrality classes. Error bars reflect the statistical and type A systematic uncertainty, the boxes represent 
the type B and C systematic uncertainties. The cancellation of uncertainties (energy scale, material budget) in the 
double ratio Ry is taken into account in the shown systematic uncertainties. 


the number of binary nucleon-nucleon collisions (A^coii = 1210.8zP 132.5, 438.4 =P42, and 77.2 zP 18 for 
the 0-20%, 20-40%, and 40-80% class, respectively) was calculated with a Monte Carlo Glauber code 


15811 using an inelastic nucleon-nucleon cross section of = 64 zP 5 mb. The decay photon spectra 


/decay Were Calculated as the product of the (/decay/?t^)|MC ratio from the decay photon calculation and 
the combined PHOS and PCM tt** spectra. Three different direct photon calculations are shown, two 
based on JETPHOX (with different parton distribution functions) 15411 . and one from Refs. oil. l52ll . 
The band around the latter reflects the factorization, renormalization, and fragmentation scale uncer¬ 
tainty whereas the bands around the JETPHOX calculations also include the uncertainty of the parton 
distribution functions. In all three centrality classes, the excess agrees with the calculated prompt di¬ 
rect photon contributions at high Pt^5 GeV/c. The contribution of prompt direct photons cannot be 
calculated straightforwardly for pj <2 GeV/c; their contribution relative to the decay photons, how¬ 
ever, is expected to be small. The excess of about 10-15% for the 0-20% centrality class in the range 
0-9 ^ Pt ^ 2.1 GeV/c indicates the presence of another source of direct photons in central collisions. 
The significance of the excess at each data point in this pT range in the 0-20% centrality class is about 
2a. Considering all data points in 0.9 < px ^ 2.1 GeV/c, the significance of fhe direef photon ex¬ 
cess is abouf 2.6a which is only slighfly larger fhan fhe significance of fhe individual poinfs due to fhe 
correlation of systematic uncerfainfies in px. 
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Fig. 4: (Color online) Combined PCM and PHOS double ratio Ry in the 0-20%, 20-40%, and 40-80% centtality 
classes compared with pQCD calculations for nucleon-nucleon collisions scaled by the number of binary collisions 
for the corresponding Pb-Pb centrality class. The dark blue curve is a calculation from Refs. |5J., which uses 
the GRV photon fragmentation function lf53l] . The JETPHOX calculations ifs^l were performed with two different 


parton disttibution functions, CTIO ll55l] and EPS09 05611 . and the BEG II fragmentation function 1157 1. 


The resulting direct photon spectra are shown in Fig. [5l Arrows represent 90% upper confidence limits. 
The same NLO pQCD calculations that were used in Fig. |4] are directly compared with the measured 
direct-photon spectra. In addition, the pQCD calculation used in the Pb-Pb direct photon prediction 
by Paquet et al. is shown as a dashed line in Fig. |51 This calculation was performed down to 
Pt ~ 1 GeV/c by using large scales p {> 2pj) and rescaling the result so that it agrees with a calculation 
done with smaller scales at higher pj. The systematic uncertainty of this calculation is estimated to be 
about 25% for Pt > 5 GeV/c, growing to about 60% at px ~ 1 GeV/c. All calculations were scaled 
with the corresponding number of nucleon-nucleon collisions in the centrality class. Similar to Ry, an 
agreement with these theoretical estimates of pQCD photon production in peripheral, mid-central, and 
central collisions for px ^ 5 GeV/c is found. An agreement between Acoii-scaled pQCD calculation and 
data for isolated direct photon yields was also found at higher px (> 20 GeV/c) by ATLAS 12711 and 
CMS Q. 


In mid-central and more clearly in central collisions an excess of direct photons at low px ^ 4 GeV /c 
with respect to the pQCD photon predictions is observed, which might be related to the production 
of thermal photons. In models in which thermal photon production in the early phase dominates, the 
inverse slope parameter reflects an effective temperature averaged over the different temperatures during 
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Fig. 5: (Color online) Direct photon spectra in Pb-Pb collisions at = 2.76 TeV for the 0-20% (scaled by a 
factor 100), the 20-40% (scaled by a factor 10) and 40-80% centrality classes compared to NLO pQCD predictions 
for the direct photon yield in pp collisions at the same energy, scaled by the number of binary nucleon collisions 
for each centrality class. 


the space-time evolution of the medium. In order to extract the slope parameter, a pj region is selected 
where the contribution of prompt direct photons is small. The pQCD contribution from the calculation 
by Paquet et al. 1^ . shown as a dashed line in Fig. [5l is subtracted and the remaining excess yield 
is fit with an exponential function oc exp(—px/Teff). The extracted inverse slope parameter is Ten = 
(297± 12®*®^±4Dy®^)MeV in the range 0.9 < px < 2.1 GeV/c for the 0-20% class and Teff = (410± 
gqstat j 4 Qsyst^ MeV in the range 1.1 < px < 2.1 GeV/c for the 20-40% class. Alternatively, to estimate 
the sensitivity to the pQCD photon contribution, the slope was extracted without the subtraction of pQCD 
photons. This yields inverse slopes of = (304 ± ID*®'±40'^^'^*)MeV for the 0-20% class and 

= (407 zh 6D*®* ± 96*^®*) MeV for the 20-40% class. The dominant contribution to the systematic 
uncertainty of the inverse slopes is due to the type B uncertainties. 


A significant contribution of blueshifted photons from the late stages of the collision evolution with high 
radial flow velocities has to be taken into account 122116311 . This makes the relation between the medium 
temperature and the inverse slope parameter less direct and a comparison to full direct photon calcu¬ 
lations including the photons emitted during the QGP and hadron gas phase is necessary to extract the 
initial temperature. A comparison to state-of-the-art direct photon calculations is shown in Fig. 0 All 
shown models assume the formation of a QGP. The hydrodynamic models, which fold the space-time 
evolution with photon production rates, use QGP rates from Ref. and equations of state from lattice 
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Fig. 6: (Color online) Comparison of model calculations from Refs. |5^62| with the direct photon spectra in 


Pb-Pb collisions at 


= 2.76 TeV for the 0-20% (scaled by a factor 100), the 20-40% (scaled by a factor 


10) and 40-80% centrality classes. All models include a contribution from pQCD photons. For the 0-20% and 
20-40% classes the fit with an exponential function is shown in addition. 


QCD. All models include the contribution from pQCD photons, however, different parameterizations are 
used. The model of van Hees et al. is based on ideal hydrodynamics with initial flow (prior to ther- 
malization) 1651] . The photon production rates in the hadronic phase are based on a massive Yang-Mills 
description of gas of n, K, p, K*, and a\ mesons, along with additional production channels (including 
anti-/baryons) evaluated with the in-medium p spectral function 11911 . Bremsstrahlung from 7i-n and K- 
K is also included 11661] . in the calculation shown here together with n-p-co channels recently described 
in Ref. i67n . The space-time evolution starts at Tq = 0.2 fm/c with temperatures Tq = 682, 641, 461 MeV 
for the 0-20%, 20-40%, and 40-80% classes, respectively, at the center of the fireball. The calculation 


by Chatterjee et al. |6lL I68l] is based on an event-by-event (2-i-lD) longitudinally boost invariant ideal 


hydrodynamic model with fluctuating initial conditions. An earlier prediction with smooth initial con¬ 
ditions was presented in Ref. 1691] . Hadron gas rates are taken from the massive Yang-Mills approach 
of Ref. 119]. Bremsstrahlung from hadron scattering is not included. The hydrodynamic evolution in 
the model of Chatterjee et al. starts at Tq = 0.14 fm/c with an average temperature at the center of the 
fireball of Tq « 740 MeV for the 0-20% class and Tq Ri 680 MeV for the 20^0% class. The calculation 
by Pag net et al. uses event-by-event (2-tlD) longitudinally boost invariant viscous hydrodynamics 
17011 with IP-Glasma initial conditions 17111 . Viscous corrections were applied to the photon production 
rates 1591 . 1721 17311. The same hadron gas rates as described above for the calculation by van Hees et al. 
are used. The hydrodynamic evolution starts at Tq = 0.4 fm/c with an initial temperature (averaged over 
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all volume elements with T > 145 MeV) of Tq = 385 MeV for the 0-20% class and Tq = 350 MeV for 
the 20-40% class. The PHSD model prediction by Linnyk et al. 16211 is based on an off-shell transport 
approach in which the full evolution of the collision is described microscopically. Bremsstrahlung from 
the scattering of hadrons is a significant photon source in this model. The comparison of the measured 
direct-photon spectra to the calculations in Fig.^indicates that the systematic uncertainties do not allow 
us to discriminate between the models. 


5 Conclusions 

The pt differential invariant yield of direct photons has been measured for the first time in Pb-Pb colli¬ 
sions at = 2.76 TeV for transverse momenta 0.9 < pj < 14 GeV /c and for three centrality classes: 
0-20%, 20-40%, and 40-80%. Two independent and consistent measurements (PCM, PHOS) have been 
averaged to obtain the final results. In all centrality classes, the spectra at high transverse momentum 
Pt > 5 GeV/c follow the expectation from pQCD calculations of the direct photon yield in pp collisions 
at the same energy, scaled by the number of binary nucleon collisions. Within the sensitivity of the cur¬ 
rent measurement, no evidence for medium influence on direct photon production at high pj is observed. 
In the low pj region, Pt ^ 2 GeV/c, no direct photon signal can be extracted in peripheral collisions, 
but in mid-central and central collisions an excess above the prompt photon contributions is observed. 
An inverse slope parameter of Teff = (297 ± 12®'^^ ±4MeV is obtained for the 0-20% most central 
collisions from an exponential function fit to the direct photon spectrum, after subtraction of the pQCD 
contribution, in the range 0.9 < pt < 2.1 GeV/c. Models which assume the formation of a QGP were 
found to agree with the measurements within uncertainties. 
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